Non-localised clustering: A new concept in nuclear cluster structure physics 
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It is shown that an angular momentum projected version of the THSR wave function, successful in 
its original form for the description of, e.g., the famous Hoyle state, has astounding properties also 
for the explanation of more compact cluster states like 20 Ne in its ground state and corresponding 
rotational bands. It is shown that these single angular momentum projected THSR wave functions 
can accurately describe the famous inversion doublet bands in 20 Ne, especially the K n — 0~ band. 
For instance, they have 99.98 % and 99.87 % squared overlaps, for 1~ and 3~ states, respectively, 
with the corresponding exact solutions of a + 16 RGM (Resonating Group Method). This finding 
sheds a completely new light on the physics of low energy nuclear cluster states in general: the 
clusters are non-localised filling the whole nuclear volume besides the one of mutual overlap where 
their probability of presence is strongly reduced. 

PACS numbers: 21.60.Gx, 27.30.+t 
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The formation of clusters is one of the most important features in light nuclei. Strongly correlated nucleons 
can compose a cluster subunit, typically, the alpha cluster. Then the relative motion between clusters becomes a 
fundamental degree of freedom, as is the single-nucleon motion in the shell model [1| . This relative motion of clusters 
displays a complex feature due to the character of the nuclear interaction and the action of the Pauli principle, i.e. 
antisymmetrisation. Refs. [2|-|5( contain recent reviews on nuclear cluster physics. 

The concept of localised clustering is a traditional understanding for the cluster structure in nuclei, which has a long 
history In 1937, Wefelmeier [6] proposed a cluster model in which the a nuclei could be considered as structure-less 
rigid alpha particles undergoing localised motion. Subsequently, Brink developed this idea and proposed a microscopic 
cluster wave function [Tj for the cluster structure in nuclei. The geometrical cluster structure was obtained by the 
energy variational method without any priori assumptions [7H2(. However, only the superposition of Brink wave 
functions, that is the GCM (generator coordinate method) Brink wave function, can lead to the exact RGM solution 
of the cluster system rather than a single Brink wave function [2tJ| . This indicates already that a single Brink wave 
function which is characterized by the localised cluster correlation cannot describe the cluster structure sufficiently 
well. 

Recently, the proposed THSR wave function lCj, |ll|, based on the concept of non-localised clustering, has brought 



a new perspective to the Hoyle state. The Hoyle state is now considered to be an alpha condensate state, in which 
the three alpha clusters occupy the same (OS) orbit and make a non- localised motion [l2j. On the other hand, also 
the 3 a RGM/GCM wave function of the ground state of 12 C was found to have about 93% squared overlap with a 



single THSR wave function 11]. Therefore, there arises the following question: Is this non- localised clustering the 



essential feature of general cluster states in nuclei? Recently, we demonstrated by using a deformed THSR type of 



wave function [13( that the concept of non-localised clustering can be extended to the more compact cluster structures 
of the positive parity ground-state band in 20 Ne [14J. The technique used angular momentum projection, already 
successfully applied in the cases of non-zero angular momentum states in 8 Be and 12 C [13J . However, in the case of 
20 Ne, where there exist negative-parity cluster states, the latter could not be included in that framework. It is, thus, 
necessary to find a new cluster wave function, inheriting the THSR spirit, for describing completely general cluster 
states in nuclei with positive and negative parity. This was the initial objective of the present work. On our way 
to the solution of this problem it turned out that we very nicely could answer in this way in addition the old and 
important question whether the cluster structure of lighter nuclei involves a more localised, crystalline picture, or 
whether the clusters are non-localised in a more fluid or gas like way. Our study will lead us to the unambiguous 
conclusion that the latter case is realised in nuclear cluster physics. 

For our purpose, we first propose a hybrid wave function which includes the Brink wave function and the THSR 
wave function as special cases. As an example, we apply this new wave function to the study of the inversion doublet 
bands in 20 Ne, especially the K v = 0j~ band. 

The hybrid wave function shall be written in the following form, 
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where 

B 2 = 2b 2 + A t (3 2 . (3) 

Here Xi, Ai, and 4>{Ci) are the center-of-mass coordinate, the mass number, and the internal wave function of the 
cluster d, respectively. ^ duster (-^i' •-•■> Rn) is the Brink wave function [7]. 

Compared with the previous THSR wave function [14j . we here introduce another generator coordinate Si in this 
new function Eq. ([1]). In this simple way, we can find that this hybrid wave function Eq. ([2]) combines the important 
traits of the Brink model and the THSR model. When S=0, Eq. (J2J corresponds to the THSR wave function and 
j3 becomes the size parameter of the nucleus. When j3 = 0, this equation is nothing more but the original Brink 
wave function and S represents the inter-cluster distance in nuclei. It should be noted that the calculation is difficult 
directly using the form of the Eq. ^ . The equivalent Eq. ([!]) adopts an integral form of the Brink wave function 



and the usual techniques for the matrix elements of the Slater determinants can be utilized, which is very useful for 
our calculations. 

Now, based on the above Eq. ([1} and Eq. (14) in Ref. 14j, the following cluster wave function of 20 Ne can be 
obtained, 



iny2 

^ No (/?,S , )ocexp( -^) 

r-S) 2 



x^4[exp(- 



5B 2 



)0(«)<X lb o)] 



(4) 



Here B 2 = b 2 + 2(3 2 , r = X\ — X 2 , Xq = (4Xi + 16X 2 )/20. X\ and X 2 represent the center-of-mass coordinates 
of the a and 16 clusters, respectively. All calculations here are performed with restriction to axially symmetric 
deformation, that is, S = (0,0, S z ). Spin and parity eigenfunctions are obtained by angular momentum projection 
technique, see below and Ref. |14J |. 
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FIG. 1. Contour map of the energy surface of the J w — + state in the two-parameter space, Sz and /3 X — Py = /3z 



TABLE I. E^in k (R) are the minimum energies at the corresponding a— 16 distance R in the Brink cluster model. -E M ^ n (/3 X , /3 Z ) 
are the minimum energies at the corresponding values of f3 x = /3 V and /3 Z in the hybrid model. E G ^ M are the GCM energies 
obtained by the hybrid model. We also show the squared overlaps between our single normalized hybrid wave functions 9^ 
corresponding to the minimum energies and the normalized GCM Brink wave functions. The squared overlaps between $ M ^ n 
and the single normalized Brink wave functions corresponding to their minimum energies are also listed. For the resonant state 
J n — 5~, we don't list the related GCM results. Units of energies are MeV. 
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We now will calculate the energy as a function of the two parameters (3 and S z . In this work, we adopt the same 
parameters as used in the previous work 14|, that is, 6 is fixed at 1.46 fm and the nuclear interaction is the Volkov 
No.l force 15]. Figure Q] shows the contour map of the energy surface of the J 71 " = + state in the two-parameter 



space, S z and j3 x — /3 y — f3 z . The minimum energy, —159.66 MeV, appears at S z = and j3 x — j3 y = /3 Z = 1.8 fm. For 
the other positive-parity states of the ground-state band in 20 Ne, the minimum-energy points also appear at S z = 0. 
Thus the hybrid wave functions are equal to the THSR wave functions in the description of these positive-parity 
states. 
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FIG. 2. Contour map of the energy surface of the J n = 1 state in the two-parameter space, S z and /3 X — ft y = j3 z 



-150- 



-152- 



a -154 



S -156 

c 

LU 

-158 ■ 



\ 1 

\ V 

\ v 

\ * '/ 
\ N ' / 
\ v ' / 

\ x / L-' 

\ x * ^ i 


/ 
/ 
/ 
/ / 

/ • \ 






+ [P=0] 

- - - r [p=o] 

+ [P=1.8] 

-— I" [3=2.41 





3 4 

SJfm] 



FIG. 3. (Color online). Energy curves of J n = + , 1 states with different widths of Gaussian relative wave functions in the 
hybrid model. 



Figure [2] shows the contour map of the energy surface of the J T = 1~ state in the two-parameter space, S z and 
fix = fly = flz- We can see that the minimum energy, —155.33 MeV appears at (3 X = (3 y = j3 z — 2.4 fm and S z = 0. 
For the J n = 3~,5~ states, the minimum-energy points also appear at S z — with different /3. As we mentioned, 
the hybrid wave function without angular momentum projection is equal to the THSR wave function at the limit 
of S z = 0, which is of positive parity. However, the hybrid wave function with odd angular momentum maintains 
its negative parity even in the limit of S z = 0, if it is normalized. This point can be explained as follows : The 
angular-momentum projection of \&Nc(/3, S)/ exp(— lOX^/b 2 ) gives us 
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FIG. 4. The energy levels of the inversion doublet bands in Ne reproduced by the hybrid wave function compared with the 
experimental levels. 



In Eq.(5), the wave function is proportional to the parameter S z . After normalisation, this parameter drops out and 
the limit S z = can be taken analytically. We, thus, obtain via the detour over the hybrid ansatz, $° M as the THSR 
wave function having definite spin L and parity independent of the localisation parameter S z . 

Of course, in practical calculations it is much more convenient for the reason given above, to take for the parameter 
S z just a very small value close to zero just within the limits of numerical accuracy. We, therefore, see that the 
introduction of the parameter S z is just a very convenient way to introduce an angular momentum projected form of 
the THSR wave function of even and odd parity. 

On the other hand, this parameter now allows us to discuss more deeply the important issue of localisation or 
non-localisation of the clusters. The competition between the two parameters /3 and S z leads to S z = 0. This is 
relevant to clarify the concept of non-localised clustering in the typical case of 20 Ne. Figure |3] shows the energy curves 
of the J 7 ' = + , l - states with different values for the widths of the Gaussian wave functions in the hybrid model. If 
/3 is fixed at 0, then the hybrid wave function becomes the Brink wave function. In this case, S z is the inter-cluster 
distance parameter. For the ground state of 20 Ne, the minimum energy appears at S z = 3.0 fm. For the J 7 ' = 1 _ 
state, the optimum position appears at S z = 3.9 fm. The non-zero values of S z seem to indicate that the a+ 16 
structure of 20 Ne favours localised clustering. This is just the traditional concept of localised clustering. However, 
now we think this argument is misleading. The non-zero minimum point S z simply occurs because the width of the 
Gaussian wave function of the relative motion in the Brink model is fixed to a narrow wave packet. If /3 is adopted 
as 1.8 fm and 2.4 fm for J 7 ' = + , l - respectively according to the minimum positions in Fig. ([1]) and ([2]), namely, 
if we use a broad enough width of Gaussian wave packet for describing the relative motion, then, we find that the 
minimum points appear at S z = in Fig. [3] It means there appears non-localised clustering. 

In Ref. 14| , we have demonstrated for the ground-state band in 20 Ne that the THSR wave functions at the minimum- 



energy points are almost 100% equivalent to the superposed Brink wave functions obtained by GCM calculations. 
The most prominent result in this Letter is that also for the K 77 — 0\ band states the hybrid wave functions at 
the minimum-energy points are almost 100% equivalent to the superposed Brink wave functions obtained by GCM 
calculations. We give in Table|T]the squared overlap values between them, which are 99.98% and 99.87% for 1~ and 
3~ states, respectively. 

Figure |4] shows the energy levels of the K 7 ' = 07/ band reproduced by the hybrid wave function, together with the 
ground-state band K n = 0^ in 20 Ne. Since it is well known that the calculation with GCM Brink wave function 
reproduced the two bands with K 77 = 1 , it is then natural that we have good reproduction of experiments also 
with our single angular momentum projected THSR wave functions. It should be noted that the results are obtained 
without any adjustable parameter. So, the good agreement with the experimental values means that this THSR-type 
wave function grasps very well the character of the relative motion of the system. 

But how shall we understand the rotational band based on the concept of non-localised clustering? In the Brink 
model, because the a- 16 Brink wave functions are deformed and parity- violating wave functions, it is easy to 
understand the reason why they can describe the inversion doublet rotational bands with K 77 = Of. On the other 
hand, in the case of the THSR wave function with S z = 0, i.e. with good symmetry, that is a wave function in the 



laboratory frame, the a and 16 clusters undergo a non-localised relative motion characterised by the large parameter 
/3 and, hence, it seems that the angular momentum projected wave function, i.e. with S z = 0, are not directly 
related with the K* — Of inversion doublet bands. In fact, the strongly anisotropic /3 X — fl y ^ (3 Z values shown 



in Table U indicate oblate shape. On the other hand, in Ref. 14j, we showed that for the ground-state band, the 
Px = Py 7^ fiz values giving the minimum energies indicate the prolately deformed shape. However, for example, 
the oblately deformed hybrid THSR wave function with /3 X — /3 y — 3.7 fm, /3 Z — 1.4 fm, and S z = fm giving the 
minimum energy for the 1~ state, has 99.98% squared overlap with the 1~ wave function projected from the prolately 
deformed THSR wave function with j3 x = f3 y = 0.1 fm, (3 Z = 3.2 fm, and S z = fm. This means that it is then 
possible to consider the prolately deformed THSR wave function as the intrinsic wave function, which can generate 
the K* — Of inversion doublet bands. A deeper discussion concerning the intrinsic state via the THSR-type wave 
function will be presented in a more elaborate version of this work. 

In summary, we treated the inversion doublet bands in 20 Ne by proposing a hybrid Brink-THSR wave function. 
From minimisation of the energy we, however, found that the spin and parity projected hybrid wave function reduces 
to a pure spin and parity projected THSR wave function, i.e. the one with S z = 0, which indicates that this separation 
parameter S z only serves in an intermediate step to produce negative parity states. We further found the surprising 
fact that these single THSR wave functions are nearly 100% equivalent to the exact RGM solution of the a+ 16 cluster 
system, which has been utilized for a long time to discuss the inversion doublet bands in 20 Ne. These results mean 
that the concept of non-localised clustering proposed by the THSR-type wave function is essential in understanding 
the cluster structures in nuclei, which is much different from the picture of separating the clusters by the localisation 
parameter in the Brink model. These astounding properties revealed by the THSR-type wave function lead to a new 
understanding of nuclear cluster physics highlighting its salient features. 
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